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Abstract 

We show that the minimal Gaussian model of nonlocal vacuum quark and quark-gluon conden- 
sates in QCD generates the non-transversity of vector current correlators. We suggest the improved 
Gaussian model of the nonperturbative QCD vacuum, which respects QCD equations of motion 
and minimizes the revealed gauge-invariance breakdown. We obtain the refined values of pion 
distribution amplitude (DA) conformal moments (£ 2N ) 7T (iV = 1,..,5) using the improved QCD 
vacuum model, including the inverse moment {x~ 1 )^, being inaccessible if one uses the standard 
QCD sum rules. We construct the allowed region for Gegenbauer coefficients 02 and 04 of the pion 
DA for two values of the QCD vacuum nonlocality parameter, A 2 = 0.4 and 0.5 GeV 2 . 
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I. INTRODUCTION 



In order to analyze meson distribution amplitudes (DAs) and form factors the generaliza- 
tion of the standard QCD Sum Rules (SRs) approach [1] has been suggested in [2, 3, 4, 5]. 
This generalization is based on the notion of the nonlocal vacuum condensates (NLC) [6, 7, 8] 
of quark and gluon fields in the nonperturbative QCD vacuum. The effects of QCD vacuum 
nonlocality appears to be very important in the pion DA analysis [9, 10, 11, 12]. 

In this approach we introduce the following gauge-invariant quark- ant iquark NLCs 1 



M s (x) = $(0)£(0,x)i>(x)) = (^) J e ax /4 da; 

o 

oo 

$(0) 7m £(0, x)jP(x)) = -ix„A [f v (a)e ax2/4 da: 
Vexp 



8(0, x) 



ig / A^dr" 
o 



:i-2) 

1.3) 



which are parameterized in the general case by distribution functions in virtualities fs(ct) 
and fv(ot), with Aq = 2a s n (ipip) 2 / 81. Explicit forms of these functions should be taken, 
generally speaking, from some concrete model of the nonperturbative QCD vacuum. This can 
be the exact solution of QCD, or some approximation, obtained, for example, in lattice QCD 
simulation. In the absence of such a model we use the first non-trivial approximation, taking 
into account only the finite value of quark momentum distribution in the QCD vacuum: 
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f v {a) = 5'( a 
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FIG. 1: Quark condensate nonlocality from lattice QCD data of Pisa group [13, 14] (solid line). 
Dashed line displays local limit, when quark condensate is constant and does not depend on distance 
z between quarks. 



1 We use the Euclidean interval x 2 = x% = —x\ — x 2 < and the subscript E will be omitted below for 
simplicity. 
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In this model, so-called "delta- Ansatz" , the average virtuality of quarks in the vacuum is 
the single parameter: 

«- <^> ' 

We have the following normalization conditions: 

poo poo ^2 

/ f s (a)da = l; / a f s (a) da = -± . (1.5) 

Higher moments of distribution /5(a) are related with higher dimensional vacuum expecta- 
tion values (VEVs) of quark fields. Delta- Ansatz (1.4) generates Gaussian form of NLC in 
coordinate representation, 

M s {x) = (^)e x > 2 / 8 ; M„(x) = ~ x„x 2 A e^ 8 . (1.6) 

For this reason below we name it as Gaussian model of NLC. The space width of this 
distribution is approximately equal to 2.5/ X q and is in a good agreement with lattice data 
(in Fig. 1 the abscissa of the symbol K corresponds to this value). This model takes into 
account one but very important property of the nonperturbative QCD vacuum — quarks 
can flow through the vacuum with nonzero momentum k and the average quark virtuality 
(k 2 ) = X 2 /2, see (1.5). It is worth to note here that Gaussian asymptotics at large values 
of \x\ differs from the anticipated exponential behavior of NLC, ~ exp(— A|x|). However, 
for the moment QCD SRs, which use averaged with the help of NLC distributions f(a) 
quantities — DA moments [2, 3], form factors [4, 15] — this wrong asymptotics of NLCs, as 
well as the more detailed information about NLC distributions, is not so important (more 
detailed discussion of this point see in [14]). 

In first papers on NLC SRs [3, 4] it was assumed that nonlocality parameters of different 
condensates (X q , Ay and XqAq) may differ. In order to simplify the NLC model and to diminish 
the number of parameters it was suggested and used in subsequent papers [10, 11, 12] to 
imply Gaussian model with the single nonlocality parameter — one and the same for scalar 
and vector NLCs (see (1.4)), and also for quark-gluon-quark (three-local) NLCs: Ay = 
XqAq — X q . As we will show in this paper, such a simplification generates the breakdown of 
transversal character of vector current correlator U^ v (q) and also of Dirac equation for the 
vector condensate (1.2). By this reason the construction of a Gaussian NLC model, which is 
minimally consistent with QCD equations of motion and minimizes the revealed breakdown 
of gauge invariance, seems to be quite reasonable. We note here that the complete restoration 
of vector current correlator transversity appeared to be impossible, because it demands to 
go outside the frames of Gaussian approximation. 

The paper has the following structure. In the next section we discuss NLCs in QCD: 
bi-local ({ip(0)^(x)) and (-^(0)7^(0;))), three-local ((ijj(0)('y5) r y f j j A u (y)ip(x))) and four-quark 
ones ((ijj(0)ip(y)'ijj(z)ijj(x))). Here we obtain an equation relating bi-local vector NLC with the 
sum of tree-locals and following from the QCD Dirac equation for the quark field operator. 
Operator product expansion for the V^V-correlator H^, taking into account nonlocalities of 
NLCs, is constructed in the third section. In the next section we analyze possible delta- 
Ansatze and find the best one, called improved Gaussian model, which minimizes the non- 
transversal part of the correlator, Hl- We show in the fifth section the results of NLC QCD 
SR analysis for the pion DA with using the improved Gaussian model. The last section 
summarizes our conclusions. 
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II. BASIC VACUUM CONDENSATES 



We use, as usual in QCD SR approach, the fixed-point gauge 



x^A^x) = 0. 



In this gauge, the gluon field operator can be expressed in terms of field-strength operators 
as follows [16] 



A;(x) = x» tci^rx) r dr . 
Jo 



A^z) dz» 



'0 

For this reason all Fock-Schwinger strings 

£(0,x) = Vexp 
if the integration path is a straight line going from to x. 



A. Bilocal quark condensates 

The vacuum expectation value (VEV) of a bilocal quark field operator can be written in 
the general form 

oo 

(nm b B(x)) = wJ{^ L /*(«) - ^A Ma)} e^da, (2.1) 
o 

where A = 2a s 7r(^ , 0) 2 /81, and functions fs(a) /y(a) parameterize the scalar and vector 
condensates, respectively. The transition to the local case is evident 

f%°(a) = 5(a) ; $ c (a) = 5'(a). 

B. Trilocal quark-gluon condensates 

It is convenient to term the quark-gluon-antiquark condensate in the fixed-point gauge 
by introducing three scalar functions Mi A 3(i 2 ,|/V 2 ) [3, 4, 11]: 

M^(x,y) = $(0) lfl (-gMy)yiP(x)) = 

= iy^Xu - ^(xy))Mi(x 2 , y 2 , (x - y) 2 ) 

+ (y l ty l/ -g^y 2 )M 2 (x\y 2 ,(x-y) 2 ); (2.2) 

M 5fMU (x,y) = (-0(O)757 M (-#A,(y))^(z)) = ie^ yx M 3 (x 2 ,y 2 , (x - y) 2 ) , 

where Ai = {— |, 2, §}A) and M 12 ,3(x 2 , y 2 , z 2 ) can be parameterized as 

oo oo 

M,(x 2 , y 2 , (x - y) 2 ) = A, [[[da, da 2 da^{a u « 2 , « 3 ) ^ WW-,) 2 )/4 _ 



000 
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We prefer to use here the hypothesis that quark and anti-quark are interchangeable in qGq- 
condensate, which means that 



fi{a u a 2 , a 3 ) = fi{a x , a 3 , a 2 ) . 
Transition to the local case for these functions is defined as follows 

f}"(<r,p,T)=6(v)6(p)6(T). 



(2.3) 



C. QCD equation of motion and nonlocal condensates 

Dirac equation for the quark field operator in massless QCD 

Viff(x) = 

allows us to write down immediately the equation of motion for the splitted quark current 

= ^(0)7 M V(a?) 

V%(x) = 0, 

where V^ 8 is covariant derivative. If we sandwich this operator equation between physical 
QCD vacuum states than we obtain the equation for condensates: 



3^0$(0) 7 ^(z)|0) = i(0|^(0) 7 ^>(x)^(x)|0); 
dPM„{x) = -iM/(x,x). 



(2.4a) 
(2.4b) 



Let us first consider the left-hand side (l.h.s.) part of this relation. By substituting (2.1) 
and delta- Ansatz (1.4) with Ay/2 = A into this part we obtain 



iAn x 2 



d^M^x) = 4 



Ax 2 



3 Ax 2 /4 



The right-hand side (r.h.s) part of (2.4) can be rewritten by using (2.2): 

oo 

- iM/(x, x) = +^|^ y* «12/ a - 9/!»(a) e^ 2 / 4 da , 

o 

where we defined the averaging ((...)) as 

1 oo 

({fi))(a)= adx da 3 fi(xa,(l-x)a,a 3 ) . 



o o 



Using (2.4), we get 



((12/2-9/!))^) e^'Ua 



3 + 



Ax 2 



3 Ao; 2 /4 



;2.s) 



We see immediately that if one uses the minimal delta- Ansatz for fi and _f 2 functions 

ff n («!, a 2 , a 3 ) = Stan- xA) 6 (cn 2 - Vi A) 6 (a 3 - z t k) , (2.6) 
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then, in order to have the same exponential in both sides of Eq. (2.5), we have to set 

xi + Vi = 1. (2.7) 

But this condition is not sufficient to fulfill Eq. (2.5): the minimal Ansatz generates only the 
first term in the square brackets in the r.h.s of (2.5), namely 3 • exp (Ax 2 /4). To cure this 
deficiency, we suggest here to use the improved delta- Ansatz: 

fi («!, « 2 , a 3 ) = (1 + Xid Xi + Yid Vi + Zid Zi ) 5 («! - x t K) 5 (« 2 - y«A) S (a 3 - *A) . (2.8) 

Then, in addition to condition (2.7), we obtain the condition for coefficients X$ and Yf. 

12 (X 2 + Y 2 )-9 (Xi + Y x ) = 1 . (2.9) 

D. Four-quark condensates 

Vacuum condensates of 4-quarks operators are usually transformed to the product of two 
scalar quark condensates by means of the Hypothesis of Vacuum Dominance (HVD) 2 

(— Tr A R\ 
-^rj M s (x 2 ) M s ((z - yf) , (2.10) 

Due to decay of correlations at large distances, say, when y 2 and {z — x) 2 are much larger 
than the characteristic scale of QCD vacuum nonlocality, 1/A 2 ~ (0.3 Fm) 2 , HVD should 
work well. In the opposite case, namely, when (z — x) 2 <^ 1/A 2 or y 2 1/A 2 , we should 
have the HVD breakdown, which is related with true 4-quarks correlations. To take this 
breakdown into consideration we can add the form factor $4 (y 2 + {x — z) 2 ), accounting for 
the separation of quark pairs (0, x) and (z,y): 

{— Ty AB\ 

$(0)Aj,(y)j>(z)B,(>(x)) = f 16N2 j M s (x 2 ) M s ({z - y) 2 ) [l + $ 4 {y 2 + (x - z) 2 )] , 

where $4(x 2 ) decreases fast for x 2 3> 1/A 2 . This modification can be done, but it does not 
appear to be very important. We will consider the influence of this modification in a separate 
paper. We suppose here, that $ 4 (x 2 ) = 0. 

III. OPERATOR PRODUCT EXPANSION OF VECTOR CURRENT CORRELA- 
TOR 

Consider now the correlator 

K = * Jd*xe^(0\T [J?(0)J+(x)] |0> , (3.1) 

of two vector currents corresponding to charged p meson 

JN(Q) = J( ) 7m (-inVof u(0) ; J+(x) = u(x) lv d{x) . 



For shortness we consider operators A and B, which include also color matrixes t a and t b 
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In the first current we have the composite operator (-inVo)^. Its action on the quark field 
is defined as 



(-mV o )%(0) = d(0) lfi (-inV y )" u(y) 



N 



where n is an arbitrary light-like vector, n 2 = 0, such that nq ^ 0. 

For shortness, we will write below Ih^, in place of IT^(g). Note that the correlator II 



l N 

depends on two vectors q and n. This dependency allows to write correlator in term s of the 
following Lorentz structures 



11^ = A N g M q u + B N g^ u q 2 + C 



N 



nq 2 



q 4 + D 



N 



nq 



q z + E 



N 



n H Qu 2 



nq 



Q 



(3.2) 



and 



= n£ [q^u-g^q 2 ] +n 



N 
T 2 



i n^n v 2 g M n v + n^q u \ 2 
9^ 9 + ( ~z^T 1 ~ ... ) 1 



+ H 



2V 
T 3 



Q/4 n u 2 

Qu Qv q 

nq 



+ n 



iV 



nq" 

q^n u + n^q u 2 

— - — q 

nq 



nq 

yrN n ti- n v 4 

nq 1 



(3.3) 



Lorentz-invariant structures An, . . . , En and 11^, IT^, IT^ are connected by the simple 
algebraic relations 



n£ = a n + d n 



N 

) 1J -T 2 



En ; II 



An + B N + D N — En ; n T = Eat — Dat 



= Aat + B N + Dr 



1 l — T -djv 1- -t^jv i n^ L — + E N — A N — B N - 
Taking into account conservation of vector current J v {x) we get: 



D r 



(3.4a) 
(3.4b) 



q v K = i"<i 2ii ? + — (n£ + n&)=o 

ng 



(3.5) 



or in terms of A/v £7 



g 1 ^ = gV (A N + B N + D N ) ^ — (CW + = 



ng 



(3.6) 



We will analyze the 11^ structure, which can be obtained using the projector n^q v j{nq). 
This structure is the most important one, because it distorts just the coefficient An- And 
namely this coefficient is responsible for distribution amplitudes (DAs) of the leading twist. 




FIG. 2: Vector quark-quark (A2vIF^,, left) and quark-gluon-antiquark (AqAqRu U , right) conden 



sates contributions to the correlator n^,. 
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FIG. 3: Four-quark condensates contributions to the correlator 11^: A 4ql n^ (left) and A 4q2 II^, 
(right). 

All 0(a s (i/>^) 2 )-terms in Yl^ u are generated by bilocal vector, quark-gluon-antiquark (see 
Fig. 2), and 4-quarks condensates (see Fig. 3): 

n£ = A 2V n^ + A^n£ + A 4Ql n* + A 4Qa n^ + (m. c.) . (3.7) 

M. C. means terms due to mirror-conjugated diagrams: for example, in Fig. 2 they correspond 
to diagrams, in which NLC are inserted in the bottom line instead of the top one. 

We are interested in the quantities corresponding to the non-transversal structure 11^: 

A,C(M 2 ) ee £L- B_ q ^ M2 k *L q = / x»<p ktL (x t M 2 ) dx 
2A nq J 

o 

with k = 2V, qAq, AQi 4Q 2 . Here, in close analogy with QCD SR approach, we work with 
Borelized quantities, which are obtained after Borel transformation B_ q 2^ M 2. Using our 
parameterizations of vacuum condensates (2.1), (2.2), we obtain for these terms: 

^ L (x,M 2 ) = -M 4 x (f v [M 2 x)-M 2 xf y (M 2 x)) ; (3.8a) 

N 2 ^2A t fff , G t (A 1 -A 2 ) jV + ^Af +2 
A qAq U L {M ) = l^-^ JJJda 1 da 2 da 3 f i {a 1 ,a 2 ,a 3 ) — - ^ ^ ; (3.8b) 

1=1 000 

A U N (M 2 ) 18 feg (A) Fl + F2) AN+2 + F3 ■ (3 8c) 

A 4Qi n L (M ) — 18 (Ar + 2)2(Ar + 3 )AA2 ■ ( 3 - 8c ) 

^4 Q2 ,l(x,M 2 ) = 36 M 2 lllMl^l ; ( 3 . 8d ) 

x 

where A = A<?/M 2 , A = 1 - A, A* = ai/M 2 , Ai = 1 - A x . Explicit form of functions F h 
Gi, and Hi are given in Appendix A. The term A 4 Q 1 TI^ r (M 2 ) is written for Ansatz (1.4). 

As we stated above, the most important for DAs of the leading twist is the coefficient An 
in front of the structure q^q^ in (3.2). Having in mind further applications for meson DAs, 
we calculate the corresponding contributions (k = 2V, qAq, 4Qi and 4Q 2 ): 

M 6 ^ Ai-n^ n' J "n l/ f 

A fc n^(M 2 ) ee — — B _q2^j\ l ,j2 ^— A = / x N <p k>T (x, M 2 )dx. (3.9) 

2A nq 2 J 

o 
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We obtain the following expressions 

V 9 2V,t(^, M 2 ) = 2M 4 xf v (M 2 x); (3.10) 

oo 

<PqAq,T( x ' M<2 ) = jjjda l da 2 da 3 f i {a 1 ,a 2 ,a 3 ) a 2 , a 3 , M 2 ) ; (3.11) 

ooo 

oo 



V9 4QliT (x, M ) = 36 JJdaida 2 fs(ai) fs(a 2 )<fi(a 1 ,a2,M ); (3.12) 

oo 

^ 4Q2iT (x, M 2 ) = 0, (3.13) 



where functions ifi(cxi, a 2 , a^, M 2 ) and tp(oci, a 2 , M 2 ) are given in the explicit form in 
Appendix A. 



IV. ANALYSIS OF GAUSSIAN MODELS 



Vector current conservation in QCD claims for the transversity (with respect to the index 
v) of the sum of contributions all quark condensates: 

An* = A 2V nf + A qAq ii? + a 4Q X + A 4Q2 n£ r + (m. c.) = o . (4.i) 

Note here that since we study Gaussian model based on delta- Ansatz (1.4), (2.8), the sum 
An* can not be equal zero exactly. The reason is very simple — we insert the Gaussian 
behavior by hands. So, the only thing we can hope to realize, is to minimize |An*| by the 
special choice of the Ansatz's parameters. More precisely, we are interested in minimization 
of conformal moments A(£ 2N )l, which are used in the QCD SR analysis of meson DAs. 
Relations between moments A(£ 2 *)l and An* are considered in Appendix B. 

In order to find these values of our parameters {X{\ we introduce the following optimiza- 
tion function (A = A 2 /(2M 2 )) 

K 

<MW) = E^<I a Ol( a ;W)| 2 >>; 

N=0 

«F(A)» = 1 f>(A = 0.024 ■ j) , 

summing up "norms" of first nontrivial K moments A(^ n )^ (A; {A,;}) with n — 0, 2, ... , 2K. 
To define the corresponding "norm" of function F(A), ((F(A))), we integrate numerically 
in A G [0 — 0.45], because just this interval of A values is physically important in QCD SRs. 
The weights wn are specified by the corresponding norms in the minimal Ansatz case [11, 12]: 

= <<|A(e 2iV > L (A; {X v = l,X i = Y i = Z 1 = 0, Xi = Vl = Zl = 1})| 2 )) . 

We introduce these weights in order to normalize contributions of different moments to the 
whole sum and to make these contributions of the same order of magnitude. For this reason 
we define wn as follows: 

$ min 

wn = z ■ ■ 

^2N 
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FIG. 4: We show functions A(£ 27V ) L (A) with N = 0,4, 10 for the improved NLC model (4.3) (solid 
line) in comparison with ones, corresponding to the minimal NLC model (dashed line). 



This receipt gives us in the minimal Ansatz case equal to unity contributions of all moments 
to the optimization function $k ({X v = 1, Xi = Yi = Z\ = 0,Xi = yi = Zi = 1}) = K + 1. 
Numerically, we use K = 5 and find 

w = 1 ; w 2 = 13 ; u> 4 = 29 ; w s = 45 ; w$ = 55 ; u>i = 59 . 



Consider now the set of available parameters {Xi} in our improved model. We apply 
delta Ansatz (1.4) with one parameter X v , relating nonlocalities in vector and scalar quark 
condensates: 



A 2 



v 



X„ A„ 



(4.2) 



For the quark-gluon-quark condensate we use Ansatz (2.8) with A = X v \ 2 q /2 and apply 
condition (2.3): 



x 



1, 2, 3). 



These parameters are not independent due to Eq. (2.9), derived from the QCD equa- 
tion of motion (2.4). After taking into account all mentioned relations we have the fol- 
lowing 7 parameters: x, X 1; X 2 , X 3 , Yi, Y 3 and X v . Minimization of the function 
$5 (x, Xi, X 2 , X 3 , Yi, Y 3 , X v ) gives us the following set of parameters 



(4.3) 



X\ = 


+0.082; 


Yt 


= Z 1 = 


-2.243; 




= x 2 


= x 3 = 


x = 0.788; X V = 1M 


x 2 = 


-1.298; 


Y 2 


= z 2 = 


-0.239; 


Ui 


= V2 


= V3 = 


l-x = 0.212; 


x 3 = 


+1.775; 


Y 3 


= z 3 = 


-3.166; 




= z 2 


= z 3 = 


1 -x = 0.212. 




Minimal Ansatz 



A(^ 2 )t 



6 ■ 



40 - 



Local Limit 



Improved Ansatz 



Minimal Ansatz 



A 



. 4 



FIG. 5: We show functions A(£°)t (left panel) and A(£ 2 )x (right panel) for the improved NLC 
model (4.3) (solid line) in comparison with ones, corresponding to the minimal NLC model (dashed 
line). 
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We will consider this set as the basic parameter set of the improved Gaussian model. 

To illustrate the quality of the improved Ansatz we show in Fig. 4 plots of the functions 
A(£ 27V )l(A) with N = 0,2,5 (solid lines) in comparison with corresponding quantities for 
the minimal Ansatz (dashed lines). As is clearly seen from this comparisons, the improved 
Ansatz (4.3) strongly suppresses the absolute values of non-transverse conformal moments 
A(£ 2Ar )L, i. e. takes the vector correlator transversity into account much better. 

In Fig. 5 we also show the moments A(£ 27V )t with N = and N = 1 for the improved 
Gaussian model in comparison with results for the minimal model [12]. 



V. PION DISTRIBUTION AMPLITUDE 



The obtained QCD vacuum model allows us to calculate moments of the pion DA 
(p n (x,fi 2 ) [17] more accurately. 



(0 | d(z) 7 M 75^(0) | tt(P)> 



= i UP" [ dxe ix{zP) <pJx,iM* 
=o Jo 



(5.1) 



The results of the analysis of (C, 2N ) n i n the NLC QCD sum rules are given in the Table I. 
One can see from this table that the values of the pion DA moments in the new Gaussian 



TABLE I: Pion DA moments (^^(/i 2 ,), determined at $ = 1.35 GeV 2 . 



Model 


U (GeV) 


N = 2 


N = 4 


N = 6 


N = 8 


N = 10 


Minimal [12] 


0.137(8) 


0.266(20) 


0.115(11) 


0.060(7) 


0.036(5) 


0.025(4) 


Ansatz (4.3) 


0.140(13) 


0.290(29) 


0.128(13) 


0.067(7) 


0.040(5) 


0.025(4) 



model of QCD vacuum are systematically different from those, corresponding to the mini- 
mal model. Allowed region for the Gegenbauer coefficients a 2 and are shown in Fig. 6. 
These coefficients define the pion DA in a form of the expansion in Gegenbauer polynomials 
C^ni^x — 1), being the eigenfunctions of the 1-loop ER-BL [18, 19] evolution kernel: 

(Pn(x; /j 2 = 1.35 GeV 2 ) = 6xx 1 + a 2 C 2 3/2 (2x - 1) + a 4 C\ l2 {2x - 1) . (5.2) 

In order to test the self-consistency of our procedure of DA restoration on the basis of 
information about its first five conformal moments, we use the same technique as in [11, 12]. 
Namely, we construct the special SR for the inverse moment (x" 1 )^ and the result of its 
processing (x _1 )^ R is compared with the inverse moment obtained from representation (5.2): 

(x _1 )° A = 3 (l+a 2 + a 4 ) . (5.3) 

For the value A 2 = 0.4 GeV 2 we get following results: 

(x _1 )° A = 3.25 ± 0.20 ; (aT 1 )^ = 3.40 ± 0.34 , 

and for the value A 2 = 0.5 GeV 2 — these: 

(s -1 )J A = 3.08 ± 0.15 ; (x _1 )^ R = 3.27 ± 0.35 . 

The obtained inverse moments in both cases are in good mutual agreement. This confirms 
the self-consistency of the pion DA recovery procedure. 
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FIG. 6: Allowed values of the pion DA parameters 02 and 04 are bounded by the solid blue line. 
Region bounded by the dotted red line represents results obtained in the minimal model [12]. Left 
panel show the results for the value A 2 = 0.5 GeV 2 , right panel - for the value A 2 = 0.4 GeV 2 . All 
values are normalized at = 1.35 GeV 2 . 



VI. CONCLUSION 

Here we considered the Gaussian model of the nonlocal vacuum quark and quark-gluon 
condensates in QCD. We analyzed the Lorenz structure of the correlator 11^ (g) of two vector 
quark currents and showed that in the minimal Gaussian model of the nonperturbative QCD 
vacuum [3, 11, 12], this correlator is non-transversal and nonlocal condensates do not satisfy 
QCD equations of motion. 

To ameliorate the situation we suggested the improved Gaussian model for nonlocal vac- 
uum quark and quark-gluon condensates in QCD, Eqs. (1.4) and (2.8). This model satisfies 
QCD equations of motion for quark fields and the revealed breakdown of gauge invariance 
is minimized by the special choice of parameters, see Eqs. (4.3). 

Using this improved model of the nonlocal QCD vacuum we analyzed QCD SRs for the 
pion DA. We revealed that in the new QCD vacuum model the NLC SRs produce again 
a 2-parameter "bunch" of admissible DAs. The allowed values of this bunch parameters 
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FIG. 7: Profiles of the pion DAs corresponding to the central points of the "bunches" for the value 
of the nonlocality parameter A 2 = 0.4 GeV 2 . Panel (a): The blue solid line represents the result 
obtained in the improved Gaussian model (symbol # on the right part of Fig. 6). Panel (b): 
The red solid line represents the result obtained in the minimal Gaussian model (BMS model [12], 
symbol o on the right part of Fig. 6). For comparison we show here also the asymptotic DA (dotted 
line) and Chernyak-Zhitnitsky (CZ) DA [20] (red dashed line). 
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a2 and 04 are shown in Fig. 6. These models are in a good agreement with the results of 
independent SR for the pion DA inverse moment, (x _1 )| R . 

We emphasize here that obtained earlier in the minimal Gaussian model of QCD vacuum 
the BMS model [12], shown in Fig. 6 by symbol o, is inside the allowed region dictated by 
the improved QCD vacuum model. This testifies to the heredity of both Gaussian models, 
the minimal one and the improved one. Moreover, that also means that all the characteristic 
features of the BMS bunch are valid also for the improved bunch: one can see in Fig. 7 
that in comparison with the CZ model [20] (dashed red line, = 0.52 and 04 = at 
fi 2 = 1.35 GeV 2 ) the NLC-dictated models are much more end-point suppressed, although 
are double-humped. 

This results in completely different values of the inverse moment: (x" 1 )^ 7, = 4.56, whereas 
in our case (a;" 1 )^ = 3.24 ± 0.20. 3 
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3 Note that smaller errors in the analysis [12] are related with their estimation only by stability with respect 
to Borcl parameter M 2 variation inside the "fidelity window" of QCD SRs. In our paper we also take into 
account internal errors of the QCD SR approach and suggest that the overall error can not be smaller 
than 10%. 
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Appendix 

APPENDIX A: 0(a s (^) 2 )-CONTRIBUTIONS 

0{a s {tj)%lj) 2 )-OTdeT terms of A^IT^, are determined by four objects. They are bilocal 
vector condensate(A 2 vn^,), 3- local quark-gluon-antiquark condensate (A^IT^,), 4-quark 
condensates {A^qJ1^ v and A 4 Q 2 n^,). We consider terms for diagrams Figs. 2 and 3. Contri- 
bution of mirror-conjugate diagrams are taken into account by symmetrical consideration, 
see Appendix B. 



(nq) 

A qAq Il» = /rfxe^/^(rf(0)7,(-mVo) JV w(0)w( 2 /) 7 pA p ( 2/ ) M ( 2/ )w(x)7^(x)); (A.2) 



A 2V nJ = — w ldxe^(u(0)l,HnV ) N d(0)d(x) ll/U (^ (A.1) 



A 4Qi n^ = ^ jdxe^ldyjdz 




x(^0)7^(-mVo)^(0)n(y)7M p (y)n(y)n(a;)7^(x)J(2)7 A A A (^)rf(^)); (A.3) 

= jdxe^jdyjdz ^ ^ 

x(rf(0)(-znVo)%n^(y) 7 M p (y)n(y)n(z)7 A A A (z)n r ^(x)7^(a;)). (A.4) 
Values A fc nf (M 2 ), k = 2V, qAq, 4Qi 4Q 2 are defined in (3.8a)-(3.8d), where 

H x = N(N+l)AxAl + H 2 ^--NH 3 , 

Ai 

H 2 = -A x ((AT + 3)AxA 2 (Ai - A 3 ) + A 3 (3A 2 + 2A 1 A 1 )) , 

tf 3 = A 2 ((JV + (N + 3)A X ) A 2 A X + A 3 (3A 2 + AiAi - (iV + 3)A 1 A 2 )) , 

Gl = - N ( N + i)A 2 2 A 1 (A 1 -A 2 ) 2 + G 2 ^-G 3 N, 

G 2 = A 2 A 2 [3(iV + l)(iV + 2)Aj -(N+1)(N + 3)A 1 A 2 2 

+ N(N + 3) AxAi A 2 + (JV + 3)AxA 2 ] 
+A!A 3 (A x - A 2 ) [(N + 1) (Ai + 2) A 2 

+ (N — 1)A 1 A 1 A 2 + 3AiA 2 + 2AiA 2 ] , 

G 3 = -A 2 (A x - A 2 ) [AiAa (A x - A 2 ) 2 + 3A 2 A 3 (A x - A 2 ) 

+A 2 Ax (JVAi + (N + 3) (AiAi + A 2 (Ax + l)))] , 
Fx = (n + 1 + A) (n + 2)(n + 3) , F 2 = A - (n + 3)[(n + l)(n + 4)A + 1] , 
F 3 = (ra + 3)A-l, 

and A = A s /M 2 , A = 1 - A, A* = a^/M 2 , Ai = 1 - A x . 

For transverse components A^II^M 2 ), see (3.9), corresponding value are given by fol- 
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lowing expressions. 

».9n xB (Ax -x) ( A -r , / xAiA 2 \ A . A . T 
) ^ (" A2Al + bg Ui - (A X - x)A a ) Al(Al ~ 

a 2 , a 3 , M 2 ) = ^ - ^ 5 (x - Ax) - (l - 5 (x - A, - A 2 ) 

x(xA 3 + A 2 (A 1 + A 3 -l)) ^_ 

^2^2 9 ( x - A 9 ( A i + A 2 - x ) ; 

^ 2 («i, « 2 , « 3 , M 2 ) = - ( 1 - ) 5 (x - Ax - A 2 ) 



A 



+ s, (2 (A, - g A 3 + A 2 (A, + A 3 - 1)) g (j _ A0 g ( Ai + ^ _ t) ; 
AxA 2 



^( ai ,a 2 ,a 3 ,M 2 ) = - * ~ ^ ^ 3 ^ 3 ~ ^ g (x - Ax) 9 (Ax + A 2 - x) 



x ((Ax - x) A 3 + A 2 (Ax + A 3 - 1)) 
AfAl 

where Aj = «j/M 2 , Aj = 1 — A; and x — 1 — x. 

/<(ai, a 2 , as) / s (a) , «x + a 2 < M 2 , a x + a 3 < M 2 2a < M 2 . (1.4), (2.8) 
, . Note that result are presented for parametric functions fi(ai, a 2 , a 3 ) and /s( a ) such 
that only «x + «2 < «x + «3 < ^ 2 2a < M 2 integration domains give contribution. 
For Ansatz (1.4), (2.8) this conditions correspond to working area of QCD sum rules. 



APPENDIX B: CONFORMAL MOMENTS 



Let us consider linear combinations of moments AIT^, 



27V 



ao l = j(2x-ir\( X )dx=^- 2 ) 2N ' k ( 2 k) J 

o fc=0 o 



x 2N ~ k (p(x) dx . 



This combinations is named conformal moments. Namely this moments are analyzed in 
QCD sum rules for meson DA. 

Reflection-symmetrical diagrams are equal to calculated diagrams. If x-density for calcu- 
lated diagrams are (fo{x), then <po(l — x) are density for M. C. diagrams and 



i 

2iV,„ f„\ J / /r. 1 \2N, 



(2x - iy N (p (x) dx = j (2x - l) 2A Vo(l -x)dx. 



That is full contribution in conformal moment of fixed diagram is equal to doubled term of 
either of diagram. 



AOl = 2 J (2x-l) 2A Vo(x)dx 
o 

Denoting 

i 

Alio = J x k (po(x) dx, 
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we immediately obtain the necessary conformal moments 

2N 



A<0 L = 2£(-2)*( 2 ^An$ 

fc=o V / 
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